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ABSTRACT 

There are no known double black hole (BH-BH) or black hole-neutron star (BH-NS) systems. We 
argue that Cyg X-3 is a very likely BH-BH or BH-NS progenitor. This Galactic X-ray binary consists of 
a compact object, wind-fed by a Wolf-Rayet (WR) type companion. Based on a comprehensive analysis 
of observational data, it was recently argued that Cyg X-3 harbors a 2-4.5 M Q black hole (BH) and a 
7.5-14.2 M Q WR companion. We find that the fate of such a binary leads to the prompt (< 1 Myr) 
formation of a close BH-BH system for the high end of the allowed WR mass (Mwr 13 M Q ). For 
the low- to mid-mass range of the WR star (M W r ~ 7-10 M ) Cyg X-3 is most likely (probability 
70%) disrupted when WR ends up as a supernova. However, with smaller probability, it may form a 
wide (15%) or a close (15%) BH-NS system. The advanced LIGO/VTRGO detection rate for mergers of 
BH-BH systems from the Cyg X-3 formation channel is ~ 10 yr -1 , while it drops down to ~ 0.1 yr _1 for 
BH-NS systems. If Cyg X-3 in fact hosts a low mass black hole and massive WR star, it lends additional 
support for the existence of BH-BH/BH-NS systems. 

Subject headings: binaries: close — stars: evolution, neutron stars — gravitation 



1. INTRODUCTION 

Ever since the discovery of the Hulse- Taylor pulsar 
1913+16 (Wcisberg & Taylor 2005), a system consisting 
of a pair of neutron stars, there has been a lot of effort to 
detect other double compact objects. So far only double 
neutron star (NS-NS) systems were found (e.g., Lorimer 
2008). In particular, no BH-NS nor BH-BH systems are 
known at the moment. 

Merging compact binaries have provided the strongest 
motivation yet to build wide-band interferometric grav- 
itational wave detectors. Second generation detectors, 
the Laser Interferometer Gravitational-wave Observatory 
(LIGO) in the US (Harry, 2010), Virgo in Italy (Acernese 
et al., 2006) and KAGRA (KAGRA 2012) in Japan, have 
the potential to observe the late inspiral and merger phases 
of compact binaries with their total mass in the range 
~ [1, 10 3 ]M Q . The range of masses they could observe 
is determined by the frequency window of ~ [10, 10 4 ] Hz 
in which they operate. Radio binary pulsars lend strong 
motivation for systems at the lower end of the above mass 
range. When LIGO, Virgo and KAGRA begin their ob- 
servations later in this decade, they could well provide the 
observational evidence for BH-NS/BH-BH systems or set 
very strong constraints on the statistics of their popula- 
tions. However, we now have a plausible progenitor of a 
BH-NS/BH-BH system in our own galaxy. 

High mass X-ray binaries (HMXBs) provide a unique op- 
portunity to study various astrophysical phenomena. At 
Warsaw Observatory we have undertaken a program to 
study HMXBs in the context of the formation of double 
compact objects like BH-BH or BH-NS systems. We have 
already provided studies of two extragalactic HMXBs that 



host the most massive known BHs of stellar origin: IC10 
X-l and NGC 300 X-l (Bulik, Belczynski & Prestwich 
2011), an analysis of Cyg X-l binary that harbors the most 
massive stellar BH in our Galaxy (Belczynski, Bulik & Bai- 
lyn 2011) and analyses of several other binaries with well 
established parameters: GX 301-2, Vela X-l, XTEJ1855- 
026, 4U1907+09, Cir X-l, LMC X-l, LMC X-3, M33 X-l 
(Belczynski, Bulik & Fryer 2012b). 

In this study, we follow the recent estimate of system pa- 
rameters for another Galactic HMXB: Cyg X-3. Zdziarski, 
Mikolajcwska & Belczynski (2012) have estimated that 
Cyg X-3 consists of a low mass BH, 2-4.5 M©, and a mas- 
sive WR star, 7.5-14.2 M©. We use this estimate to cal- 
culate the future evolution of Cyg X-3 to check whether 
this binary may provide any observational constraints on 
as-yet undetected BH-BH and BH-NS systems. The past 
evolution of Cyg X-3 was studied in detail by Lommcn ct 
al. (2005). 

2. ESTIMATES 

2.1. The future evolution of Cyg X-3 

In this section we describe in detail Cyg X-3 evolution- 
ary scenarios leading to the formation of a BH-BH or 
BH-NS system. For demonstration we choose extremes 
of the allowed mass range of BH and WR components: 
M BH i = 2.0 with M W R = 7.5 M Q and M B hi = 4.5 
with Mwr = 14.2 M Q . Additionally, we evolve the most 
likely configuration Mbhi = 2.4 with Mwr = 10.3 M© as 
adopted from Zdziarski et al. (2012). First, we employ 
our standard model for stellar and binary evolution (see 
below) and in the next section we will test how various 
evolutionary uncertainties change our results. 
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We evolve three M zams = 27, 36, 50 M single stars 
(Hurley, Pols & Tout 2000) at solar metallicity Z = 0.02 
with our calibrated wind mass loss rates (Belczynski et al. 
2010a). At t — 6.5, 5.2, 4.4 Myr these stars leave the main 
sequence (MS) and become Hertzsprung gap objects with 
helium core of M\yR,i = 7.5, 10.3, 14.2 M , respectively. 
We expose this core right after MS and make them naked 
helium stars: massive WR objects. We pair these WR 
stars with M B m = 2.0,2.4,4.5 M Q BHs and create three 
binaries all with orbital period of P rb = 4.8hr as observed 
for Cyg X-3. Such an approach allows a WR star to shed a 
maximum amount of mass, reducing its chance to become 
NS or BH, and thus makes our estimates conservative. We 
should note that this simplified picture in which the WR 
star is evolved independently of its companion ignores the 
effect of past interactions with the BH companion on the 
future evolution of the WR star, which could change the 
rate of future mass loss. However, the mass loss rates that 
we apply to WR companion are derived based on a large 
population of WR stars, including those in binaries (see 
below) . 

The binary separations at the currently observed or- 
bital period of the three synthetic binaries are a = 
3.0, 3.4, 3.8 R and the Roche lobe radii of WR compo- 
nent are i?i bc = 1-5, 1.7, 1.8 R©. The corresponding radii 
of WR stars are i?wR = 0.8, 1.0, 1.2 R . Massive helium 
stars do not expand at any significant level during evolu- 
tion. Therefore, no Roche lobe mass transfer episode is 
expected in the future evolution of these systems. 

The lifetimes of these WR stars arc t\u e = 0.95, 0.77, 0.64 
Myr and at the end of their evolution the masses drop to 
MwR,f = 5.5, 6.8, 8.2 M . We have applied wind mass loss 
rates adopted from Hamann & Koesterke (1998) that take 
into account effects of clumping 
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This prescription is based on detailed stellar evolutionary 
calculations combined with comprehensive WR wind mod- 
els and calibrated on observations of three WR stars. We 
obtain WR star luminosity (L) from evolutionary models 
of Hurley et al. (2000). 

At the time of explosion, orbits have expanded to 
a = 3.9, 4.7, 5.6 R due to the WR wind mass loss 
and WR stars have formed massive CO cores (Mco = 
4.0,5.0,6.2 M ). The WR stars undergo a core collapse. 
We use Fryer et al. (2012) rapid explosion model that re- 
produces the observed mass gap between neutron stars 
and black holes (Belczynski et al. 2012a). For initial 
Mwr,) = 7.5, 10.3 M WR stars the core collapse is fol- 
lowed by Type lb supernova and neutron stars form with 
mass Mns = 1.5, 1.7 M , respectively. The mass loss and 
natal kicks associated with supernova are very likely to 
disrupt these two synthetic binaries (/disruption = 0.69). 
However, there is also a significant chance to form a close 
BH-NS system. We estimate the probability of the BH-NS 
system formation with gravitational merger time T merger 
shorter than 10 Gyr at the level / c i OS c = 0.14, 0.12 for the 
low and intermediate mass realization of Cyg X-3. We 
have applied a distribution of natal kicks as derived from 
the population of Galactic pulsars (Hobbs et al. 2005). 
The kicks have a random direction and their magnitude is 



taken from a single Maxwellian with a = 265 km s _1 . Ad- 
ditionally, we lower the magnitude of natal kicks due to the 
amount of fall back expected in each supernova explosion 
Mdck = (1 - /fb)Vkick (Fryer et al. 2012). The fall-back 
amount, the fraction of matter that was initially ejected in 
a supernova explosion but that is later accreted back onto 
the compact object, was estimated to be /fb = 0.14,0.16 
in case of models with initial Af\yR,i = 7.5, 10.3 M , re- 
spectively. 

For the high component-mass realization (Mwr,; = 
14.2 M ), within our rapid supernova explosion model the 
entire WR star falls into a BH (i.e., /fb = 1.0). We 
assume the loss of 10% of the gravitational mass dur- 
ing the collapse through neutrino emission, which leads 
to a slight orbital expansion and a small induced eccen- 
tricity (a = 6.0 R , e = 0.07). In the end the second 
BH forms with mass Mbh,2 = 7.4 M . No natal kick 
is imparted on a BH within this model. Therefore, a 
close BH-BH system forms with a chirp mass of M Cj d co = 
(MiM 2 ) 3/5 (Mi + M 2 )~ 1/5 = 5.0 M Q and a merger time 



of T„ 



0.5 Gyr. In this particular case with no natal 



kick there is no other possibility (/ c ioso = 1-0), within the 
framework of our model, then to form a close (Tmorgcr < 10 
Gyr) BH-BH system. 

If we conservatively assume that at any given time there 
is only 1 such system as Cyg X-3 it means that Galactic 
birth rate is at the level 7?-birth = 1/twR- As listed above, 
lifetimes of the considered massive WR stars are very short 
^wr Ss 1 Myr. Since the merger times of systems that we 
include in our analysis are relatively short (T mGrgcr < 10 
Gyr) and the star formation was approximately constant 
in the Galactic disk over a long period of time (~ 10 Gyr) 
the Galactic merger rate of BH-NS and BH-BH systems 
may be estimated from 
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we discuss statistical uncertainties in the estimate of the 
birth rates in section 2.4. For our three Cyg X-3 real- 
izations we obtain Rmw = 0.19,0.20 Myr" 1 in the cases 
of BH-NS formation and Kmw = 1-56 Myr -1 in the case 
of BH-BH formation. The main results presented in this 
section are summarized in Table 1 (top 3 entries) and il- 
lustrated in Figure 1 (top panel). 

2.2. Range of evolutionary uncertainties 

Alternatively to our standard approach, we employ WR 
wind mass loss rates adopted from Zdziarski et al. (2012): 

/ Af \ 2.93 
{dM/dt) = 1.9 x 10~ 5 —^7- M © yr" 1 - ( 3 ) 
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Zdziarski et al. (2012) selected data points from Nugis & 
Lamcrs (2000) for WR stars of similar mass and type as the 
one in Cyg X-3. Nugis & Lamers (2000) have estimated 
clumping corrected mass loss rates based on observations 
and modeling of 64 Galactic WR stars. We refer to this 
wind mass loss prescription as empirical as it is heavily de- 
rived from observations in contrast to the method used in 
our standard approach (mostly originating from theoreti- 
cal predictions). In particular Zdziarski et al. (2012) em- 
ployed data points from Table 5 of Nugis & Lamers (2000) 
for WN stars. It is worth noting that the majority of these 



3 



WN stars (19 out of 34) are in binary systems. The results 
of calculations with those winds are presented in Table 1 
(models marked with "Wind: empiri" entry) and in Fig- 
ure 1 (bottom panel). It is noted that both wind mass 
loss prescriptions give very similar results. There are no 
significant changes in our analysis. There is only a slight 
increase of BH mass (A/bh.2 = 8.0 M Q ) in the case of BH- 
BH formation. 

Next, we alternate our approach to natal kicks. Wc 
adopt high kicks for all compact objects. It means that 
both neutron stars and black holes receive the kicks from 
the same distribution described by a Maxwellian with 
a = 265 km s _1 . There is no kick decrease factor ap- 
plied in this case. Results for such an approach are listed 
in Table 1 (marked with "Kicks: high" entry). There is 
no significant changes for the BH-NS formation, as neu- 
tron stars in our standard model analysis were receiving 
almost full kicks (low fall back). There is a noticeable de- 
crease in the formation efficiency of close BH-BH systems 
(/close = 0.68) that leads to corresponding decrease in the 
Galactic merger rate (Kmw = 1-06 Myr -1 ). The decrease 
is caused by binary disruptions due to high natal kicks 
applied to black holes in this model. However, the for- 
mation efficiency is still significant despite the rather high 
kicks that black holes receive in this model, because the 
binary's initial orbital velocity, ~ 800 km s _1 , is already 
large compared with typical kicks. 

Our calculations so far were based on the rapid super- 
nova explosion engine, in which a sufficiently massive pro- 
genitor does not explode at all, gets no natal kick and 
forms quite massive black hole. In the rapid model this is 
the reason for the emergence of a mass gap between NS and 
BH masses (Belczynski et al. 2012a, although see Kreid- 
berg et al. 2012). On the one hand, lower mass stars are 
subject to strong explosions and form neutron stars (as in 
case of Mwr = 7.5, 10.3 M Q ). On the other hand, higher 
mass stars do not explode at all and form massive black 
holes (as for AfyvR = 14.2 M Q ). This explains the rather 
sharp transition for all models considered thus far from 
Mwh = 10.3 M Q forming a neutron star (Mns = 1.7— 
1.8 Mq) to Mwr = 14.2 M forming quite massive black 
hole (M B H2 = 7.4-8 M ). 

As a final variant, we have modified our approach to 
supcrnovae explosions. We now use the delayed explosion 
model of Fryer et al. (2012). This model generates a 
continuous compact-object mass spectrum (i.e., without a 
mass gap) . The results for this set of calculations are listed 
in Table 1 (marked with "SN: delayed" entry). Due to the 
delayed nature of the explosion engine, neutron stars are 
typically more massive as proto-neutron stars have more 
time for accretion between the core bounce and actual ex- 
plosion. The explosions are also typically less energetic, 
since at later times, after some cooling, there is less en- 
ergy to drive the explosion. This leads to two changes 
for BH-NS formation. First, there is a noticeable increase 
in the formation efficiency (lower explosion energy leads 
to increased fall back and smaller natal kicks). Second, 
we note a significant increase in the NS mass. In fact, 
for the most likely Cyg X-3 configuration (Mbhi = 2.4; 
M WR = 10.3 Mq), instead of forming a NS (M NS = 1.7- 
1.8 M Q in all other models), wc obtain a low-mass black 
hole (M B H2 = 2.7 M Q ) and thus transition from BH-NS 



to BH-BH formation. In the case of the Cyg X-3 high- 
mass configuration we predict as before the formation of 
a BH-BH binary. The black hole in this model is formed 
with moderate fall back (/fb = 0.43), and therefore re- 
ceives a moderate natal kick. It has significantly lower 
mass (Mbh2 = 3.9 M©), since the majority of the ejected 
material was not a subject to fall back. As a result, the 
formation efficiency is significantly decreased (/ c iosc = 0.5) 
relative to all other models. Figure 2 shows the mass of 
the final compact object as a function of the initial mass of 
the WR star for all alternative models considered in this 
section. 

2.3. LIGO/VIRGO detection rate estimate 

In order to compute the advanced LIGO/ Virgo event 
rates, we use the same detection threshold that was 
adopted in (Abadie et al., 2010) - a signal-to-noise ratio 
(SNR) of 8 in a single interferometer at the sensitivity of 
advanced LIGO. This is a simplified treatment, since ac- 
tual detector sensitivity will depend on the network config- 
uration, data quality, and the details of the search pipeline, 
but the uncertainties introduced by this simplifying as- 
sumption are smaller than the uncertainties in the merger 
rates. 

Previous calculations in the population-synthesis litera- 
ture (e.g., Belczynski et al., 2010b) were based on a sim- 
ple scaling of the detection volume for arbitrary systems 
with the horizon distance do for a NS-NS binary (the dis- 
tance at which an optimally oriented source with compo- 
nent masses Mi = M 2 = 1.4 M could be detected at an 
SNR of 8): 

v 4tt / d Q \ 3 f M cAco \ 15/ \ ... 

^LIGO = Pgal-- ~ f -J-. T^MW- (4) 

*J \ Jpos / \ ^^c.nsns / 

Here, p ga \ is the density of Milky Way-like galaxies, the 

chirp mass M c , nsns = M 1 3/5 Af| /s (Mi +M 2 )~ 1/5 = 1.2 M 
and the correction factor / pos = 2.26 takes into account 
the non-uniform pattern of detector sensitivity and ran- 
dom sky location and orientation of sources (Finn, 1996). 

However, this simple calculation suffered from several 
shortcomings. The sensitivity distance scales as d ~ 
^(A^cdco/A^cnsns) 5 ^ 6 only for masses sufficiently low 
that the gravitational- wave signal spans the bandwidth of 
the detector. The waveforms used to compute the horizon 
distance included only the contribution from the inspiral 
portion, not the merger and ringdown signals. Finally, cos- 
mological effects due to the expansion of the universe have 
not been included: the rcdshifting (dilation) of masses in 
the gravitational waveform; the difference between the vol- 
ume as a function of luminosity distance d and comov- 
ing volume; and the difference in the rate of clocks in the 
source and merger frames. Here, we include these effects 
by using the following procedure when integrating over 
shells centered on the detector: 

• For each shell at a given rcdshift z, compute the lu- 
minosity distance d(z) and the comoving volume of 
the shell dV c (z) using standard cosmology (see, e.g., 
Hogg, 1999) with n M = 0.272, fl A = 0.728, h = 70 .4 
(WMAP 7 results). 

• Given a particular combination of compact-object 
masses from Table 1, compute the SNR for an opti- 
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mally oriented source at distance d(z), 

where we use the inspiral-merger-ringdown waveform 
family IMRPhenomB (Ajith, 2010) to compute the 
frequency-domain waveforms h(f) using redshiftcd 
component masses M — > M(l + z) and the zero- 
detuning, high-power Advanced LIGO noise power 
spectral density S n (f). 

• Compute the fraction of detectable mergers fdet(z) 
from this shell by considering the fraction of sources 
for which the projection function 0(t, 9, ?/>, <fr) [given 
in Eqs. (3.4.b-d) of (Finn, 1996)] is large enough so 
that (<d/4)SNR > 8. The cumulative distribution 
function of O can be computed numerically (the ap- 
proximate expression given by Finn, 1996, not being 
sufficiently accurate) via a Monte Carlo over the in- 
clination angle i (whose cosine is uniform in [0, 1]), 
the sky-location spherical coordinates 9 (whose co- 
sine is uniform in [0, 1]) and <fi (uniform in [0,27r]), 
and polarization ip (uniform in [0,7r]). 

• The contribution of the shell to the detection rate is 
given by fdct(z) dV c ( z ) PgaiT^-M w ( 1 + z ) ~ 1 . We assume 
that the intrinsic merger rate in the source frame is 
independent of redshift (i.e., we don't include the 
redshift dependence of star formation rate or metal- 
licity). Thus, we use Kmw from Table 1, and scale 
it by the space density of Milky Way-like galaxies, 
for which we use 0.01 per Mpc 3 of comoving volume. 
The final factor of 1/(1 + z) reflects the time dila- 
tion between the source clock, used to measure the 
merger rate, and the clock on Earth, used to measure 
the detection rate. 

We can now integrate over multiple shells (or, in prac- 
tice, compute an approximate Riemann sum over the shells 
numerically) to obtain the advanced LIGO/Virgo detec- 
tion rates that we report in Table 1. Including merger and 
ringdown waveform phases and accounting for cosmologi- 
cal redshift corrections decreases the detected event rates 
by > 15% for BH-NS systems and by £ 30% for BH-BH 
systems relative to the simple scaling of Eq. Qj, largely be- 
cause the comoving volume within a luminosity distance 
d is significantly smaller than (4/3)-7rd 3 at non-trivial red- 
shifts. 

2.4. Statistical uncertainty on merger and detection rates 

We note that the rates derived above do not account 
for the statistical uncertainty associated with observing a 
single binary like Cygnus X-3. In practice, we do not know 
the exact rate from a single observation even if / c i OS g and 
^wr are known perfectly in Eq. [2] Assuming that the birth 
of binaries like Cygnus X-3 is a stochastic Poisson process 
with rate 7?-birth, the probability of electromagnctically ob- 
serving exactly one such system is 

p(l obs.|^birth) = ^bhthtwR exp (-T^birthtwR)- (6) 

We can compute the probability distribution on the rate 
given a single observation by using Bayes' theorem: 

p(Tl h nth\l obs.) oc p(7?. b irth)p(l obs.|7?.birth), (7) 



where p(lZbirth) is the prior probability on the birth rate 
of such systems. 

If a flat prior is chosen, p(TZbinh) = const, the most 
likely birth rate is 1/twR, the value used in Eq. How- 
ever, the 90% credible interval on the birth rate extends 
from 0.35/^wr (at the 5th percentile) to 4.74/iwR (at the 
95th percentile). 

Meanwhile, if an uninformativc Jeffreys prior on the rate 

— 1/2 

is chosen, p(1Zbhth) cx lZ hirth , the most likely birth rate is 
halved to 1/2/<wr- The 90% credible interval is shifted 
downward to 0.17/iwR ~~ 3.9/twR- 

Finally, we can consider the case where we assume 
that there is at least one rather than exactly one bi- 
nary in the Galaxy in the same stage as Cygnus X- 
3 (i.e., there may be other similar Galactic systems 
which have not been observed). In that case, we are 
interested in p(lZbh-th\> 1 detection), which scales with 
p(> 1 detection 1 7?.birth) = 1 — cxp (— 7?-birthtwR)- For the 
Jeffreys prior, the posterior birth rate distribution peaks 
at 1.26/^wr- 

The merger rate in all cases is given by TZuw = 
7^-birth //close, and so ranges by the same pre- factors rel- 
ative to the rate given in Eq. [5] The large range of purely 
statistical uncertainty, which spans about a factor of 5 
above and below the value in Eq. [21 reflects the difficulty 
of making robust inference from a single observation. 



3. DISCUSSION 

We have calculated the future evolution of Galactic bi- 
nary Cyg X-3 harboring a compact object and a WR star. 
We have employed a recent Cyg X-3 component mass es- 
timate from Zdziarski et al. (2012) and following their 
arguments we have assumed that the compact object in 
Cyg X-3 is a black hole. 

Our results indicate that the future evolution and fate 
of Cyg X-3 is a strong function of mass of the WR star. 
Within the measurement uncertainties Cyg X-3 may ei- 
ther form a close BH-BH binary at the high end of the 
allowed WR mass (Mwr ~ 14 M Q ), or form a close BH- 
NS (~ 15%), a wide BH-NS system (~ 15%) or get dis- 
rupted producing single BH and NS (<~ 70%) at the low 
end and middle of the allowed mass range for the WR star 
(M WR ~ 7-10 M ). 

We have estimated the advanced LIGO/VIRGO detec- 
tion rates in case of the close BH-BH and BH-NS for- 
mation. The rates are significant: in the case of BH-BH 
formation 7^ligo ~ 10 yr . This is the first empirical 
estimate of a BH-BH detection rate based on a Galactic 
system. Previous empirical estimates were based on ex- 
tragalactic high mass X-ray binaries in small star forming 
galaxies IC10 and NGC300 (Bulik et al. 2011). The pre- 
dicted rates that are based on observations within Milky 
Way (high mctallicity) are much lower than estimated for 
the above two low-metallicity galaxies. The BH-BH detec- 
tion rate extrapolated from IC10 X-l and NGC300 X-l for 
advanced LIGO/Virgo is TvLligo ~ 2000 yr -1 , where we 
used conservative mass estimates and mean merger rates 
from (Bulik et al. 2011) and converted merger rates to 
detection rates as discussed in section 2.3. The low metal- 
licity environment can significantly boost close BH-BH for- 
mation rates as explained by Belczynski et al. (2010b). 
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The detection rate for BH-NS systems forming via the 
Cyg X-3 channel are low, 7£ligo ~ 0.1 yr -1 . However, our 
rates are only lower limits as more binaries similar to Cyg 
X-3 may currently exist in the local Universe and close 
BH-NS systems may potentially form via other formation 
channels. Among about ~ 200 Galactic and extra-galactic 
high-mass X-ray binaries only a handful have established 
parameters (Liu, van Paradijs & van den Hcuvel 2005, 
2006) allowing for detection rate prediction (Bclczynski 
et al. 2012b). Recent population synthesis calculations 
provide rates that are typically a few detections per year 
for BH-NS systems for advanced detectors (Dominik et al. 
2012; Belczynski et al. 2012c). Our prediction is only the 
second empirical estimate for BH-NS detection rates. The 
first one was obtained for another Galactic system Cyg X- 
1 (Bclczynski et al. 2011). We note that our current rate 
(~ 1 detection per decade) is 10 times higher than the one 
obtained for Cyg X-l (~ 1 detection per century). 

Beyond our standard evolutionary calculations we have 
performed several models to check the validity of our con- 
clusions. We have varied WR wind mass loss rates, na- 
tal kicks that compact objects receive in supernovae, and 
the supernova explosion mechanism that alters the NS/BH 
mass spectrum. As long as we stay within the framework 
of rapid supernova explosion mechanism there are no sig- 
nificant changes to our conclusions. We find a range of de- 
tection rates 7?-ligo = 0.09-0.15 yr _1 for BH-NS systems 
and7?-LiGO = 7.7-12.4 yr _1 for BH-BH systems. Our stan- 
dard supernova model, which employs rapid explosions, is 



motivated by the existence of the mass gap between neu- 
tron stars and black holes (e.g., Bailyn 1998; Ozel et al. 
2010; Farr et al. 2011). However, if the mass gap is not 
an intrinsic signature of the BH/NS mass spectrum but 
is caused by some observational bias as recently claimed 
by Krcidberg et al. (2012), our conclusions change. For 
the delayed supernova model, which is consistent with the 
absence of a mass gap, we find that BH-NS formation 
occurs only at the lowest allowed mass for the WR star 
in Cyg X-3 (Mwr ~ 7 M ) with a slightly higher de- 
tection rate 7?-ligo = 0.17 yr _1 . BH-BH formation is 
found in a broader mass range allowed for the WR star 
(M WR ~ 10-14 M Q ). The second BH falls right within 
the mass gap (Mbh2 = 2.7-3.9 M©) and the rates are sig- 
nificantly smaller, 7^-ligo = 0.41-2.8 yr _1 . Finally, we 
note that all merger and detection rates have statistical 
uncertainties of approximately a factor of 5 in cither di- 
rection due to the limited observational sample of a single 
system. 
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Table 1 
The Fate of Cyg X-3 
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a Either theoretically (eq. [T]) or empirically (eq. [3]) based WR wind mass loss rates are applied. 

b High: NS and BH kicks are taken from a Maxwellian with a = 265 km s _1 . Low: the high kicks are decreased 
proportionally to the amount of fall back for both NSs and BHs. 

c Compact-object formation model: either via rapid supernovae (mass gap) or delayed supcrnovae (no mass gap). 
d Type of binary formed followed by the mass of compact objects. 

e We only count binaries that have formed with merger time shorter than 10 Gyr (others are disrupted or form wider 
systems). 
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Fig. 1. — Evolutionary prediction for Cyg X-3 WR star (MwR.i = 10.3^2 g M©). Top panel shows the wind mass 
loss rate based on theoretical calculations, while results in bottom panel are based on observationally estimated WR 
mass loss rates. Independent of the adopted wind mass loss rate the WR component becomes either a neutron star for 
low- to medium-mass progenitors (Mwr,i = 7.5-10.3 M©), and a black hole at the high end (MwR,i = 14.2 M©) of the 
allowed WR mass range. The mass of the WR component at the end of its evolution is marked with Mwr f, while the 
mass of compact object formed after core collapse/supernova explosion is marked with A/ns/bh- Compact objects with 



M- 



NS/BH 



< 2 M Q are assumed to be neutron stars, and above that black holes. 
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Fig. 2. — The dependence of the second compact object mass on the current WR star mass for Cyg X-3. There is a sharp 
transition from NS to BH formation for the rapid SN engine at around Mwr,i ~ 12-13 M Q . The delayed SN engine allows 
for a steady increase of compact-object mass and the NS/BH transition depends sensitively on the unknown maximum 
NS mass. The assumed kick velocity has no impact on the compact-object mass, so two kick models (theory/low/rapid 
and theory /high/rapid) yield the same mass spectrum for the second compact object. 



